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ABSTRACT 

XMM-Newton observations of the bright Seyfert 1.9 galaxy MCG-5-23-16 have revealed a broad FeKa 
emission line which is nearly symmetric in contrast to the broad and red-shifted asymmetric FeKa line 
sometimes observed from Seyfert 1 galaxies. The FeKa line has two distinct components - a narrow 
unresolved component with equivalent width of ~ 40 eV and a broad component with full width at half 
maximum of ~ 40000 km s -1 and equivalent width of ~ 120 eV. An absorption feature at ~ 7.1 keV 
has also been observed. The energies of the emission and absorption features are consistent with those 
arising from neutral iron. The broad component is consistent with an FeKa emission line expected from 
a relativistic accretion disk around a Schwarzschild or a Kerr black hole. Alternatively, most of the flux 
in the broad component could also be modeled as reflection emission which mimics emission line like 
feature due to the presence of iron K-shell edge at ~ 7.1 keV, however, the reflection fraction, R ~ 3, is 
much higher than that inferred from the BeppoSAX observations (R ~ 0.5). The disk inclination angle 
of ~ 47deg, inferred from the disk-line fits, and the absorption column (Nh ~ 10 22 cm -2 ), inferred from 
the low-energy spectral curvature due to photoelectric absorption, suggest that our line of sight passes 
through the outer edge of a putative torus and are consistent with those expected for a Seyfert 1.9 galaxy 
falling within the unification scheme. The strength of the narrow iron Ka emission and the optical depth 
of the iron K absorption edge suggest their origin in the putative torus with JVh ~ 10 cm -2 in the inner 
regions and JVh ~ 10 22 cm -2 in the outer edges. The strength of the broad component of FeKa varied 
by a factor of ~ 2 between the two XMM-Newton observations taken ~ 6 months apart, while the narrow 
component of FeKa and the continuum flux did not appear to vary appreciably. There is evidence for a 
weakening in the strength of the broad iron Ka emission with the flattening of the observed continuum. 
This can perhaps be explained if the shape of the continuum is coupled with the ionization stage of the 
reflector. 

Subject headings: Active galactic nuclei: accretion disks — line: formation 



1. INTRODUCTION 

The Seyfert-type active galactic nuclei (AGNs) appear 
to be intrinsically X-ray loud and the variation in their X- 
ray emission is largely due to different absorbing columns 
along the respective lines of sight. The type 1 and type 2 
Seyferts and the intermediate Seyferts (types 1.5,1.8, and 
1.9), all show a similar form in their intrinsic X-ray emis- 
sion - a power law with photon index of ~ 1.9 which can 
extend up to a fewxlOO keV, an FeKa line at ~ 6.4 keV, 
and a reflection hump in the ~ 10 — 100 keV region. How- 
ever, the obscured Seyferts show a low-energy cut-off of 
the intrinsic power law due to photoelectric absorption. 
The exact energy of the cut-off depends on the depth of 
the absorption column. The obscured Seyferts also show 
unabsorbed soft X-ray emission composed of a number 
of photo-excited emission lines (Kinkhabwala et al. 2002; 
Schurch et al. 2003). This component is often extended, 
non- variable and cannot be the emission from the accretion 
disk. X-ray emission arising from the accretion disks of ob- 
scured Seyfert galaxies (e.g., the broad iron Ka emission) 
provide us with the opportunity to study the geometry 
and physical conditions in the inner regions surrounding 
a super massive black hole (SMBH), thus helping us in 
understanding the diversity in the AGN characteristics. 

The iron Ka emission is the most prominent ubiquitous 
feature in the 2—10 keV region of the X-ray spectra of 



AGNs. This emission appears to be a general character- 
istic of AGNs. Complex iron Ka emission consisting of 
multiple components and/or asymmetric profiles arising 
from neutral and/or ionized iron have been observed from 
many Seyfert 1 galaxies and provide the strongest evidence 
for the existence of cold material in the central regions of 
AGNs. The broad iron Ka emission is the key feature in 
studies of the geometry and nature of the accretion disks 
around black holes. The line profile can also be used to in- 
fer the Schwarzschild or Kerr nature of black holes. While 
some Seyfert 1 galaxies show broad iron K-shell emission, 
there are only a few Seyfert 2 galaxies with broad iron Ka 
emission. This is likely due to the difficulty in detecting 
the broad wings in the presence of large absorbing columns 
along our line of sight to a type 2 nucleus. The narrow- 
emission line galaxies (NELGs) are good targets for study 
of the iron Ka emission from the obscured AGNs. Op- 
tically, NELGs are the same as Seyfert 2s and show nar- 
row optical emission lines, often with broad wings of Ha 
and P/3 which makes them Seyfert 1.9 nuclei. NELGs are 
bright and variable X-ray sources which were discovered 
in early X-ray surveys (Marshall et al. 1979; Griffiths et 
al. 1978). In this paper, we present a detailed study of the 
iron Ka emission from an obscured Seyfert galaxy MCG- 
5-23-16 using XMM-Newton observations. 

MCG-5-23-16 is an SO galaxy (Ferruit, Wilson, & 
Mulchaey 2000) first discovered in X-rays by Schnopper et 
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al. (1978) who also found its' optical spectrum to be that 
of a narrow emission-line galaxy at a redshift of 0.0083. 
Veron ct al. (1980) classified this galaxy as a Seyfert 1.9 
based on the possible presence of a broad component in 
Ha. MCG-5-23-16 with a V-magnitude of 13.7 is one of 
the brightest Seyfert galaxies in hard X-rays and is also one 
of the few Seyferts detected at 7-ray energies (Bignami et 
al. 1979; Pollock et al. 1981). Its 2-10 kcV flux, measured 
to be - 8 x lO" 11 erg cm" 2 s" 1 in 1978 (Tennant 1983), 
decreased by a factor of four in 1989 (Nandra & Pounds 
1994), and then increased to <~ 9 x 10~ n erg cm~ 2 s _1 
in 1996 (Weaver, Krolik, & Pier 1998) before decreasing 
again to 7 x 10~ n erg cm~ 2 s _1 in 2001 (this work). Pre- 
vious X-ray observations of MCG -5-23-16 showed a highly 
absorbed X-ray spectrum (Nn ~ 10 22 cm -2 , Turner & 
Pounds 1989) and a strong iron Ka fluorescence line with 
equivalent width of 300 - 400 eV (Singh, Rao, & Vahia 
1992; Nandra & Pounds 1994; Weaver et al. 1997; Weaver, 
Krolik, & Pier 1998). Weaver et al. (1997) reported a 
complex iron K-shell line profile in an ASCA observation 
of this source. These authors discovered a narrow core at 
the systemic velocity of the galaxy and wings to the red 
and blue sides of this core. They modeled the line profile 
with three Gaussian having central energies 5.37, 6.37, and 
6.58 keV. Alternatively, the data could also be modeled 
with a line profile from an accretion disk, superimposed 
upon a narrow Gaussian profile. These data have pro- 
vided some of the best evidence to date for contributions 
from two distinct regions to the observed iron line profile. 

In this paper we present the first results from the XMM- 
Newton observations of MCG-5-23-16. The paper is orga- 
nized as follows. In Sect. 2, we describe the XMM-Newton 
observations and the data selection. In Sect. 3, we present 
the detailed spectral analysis of the FeKa region. We dis- 
cuss our results in Sect. 4 followed by conclusions in Sect. 5. 

2. OBSERVATIONS AND DATA SELECTION 

XMM-Newton observatory (Jansen et al. 2001) has three 
Wolter type 1 X-ray telescopes with three European Pho- 
ton Imaging Cameras (EPICs) - one PN (Struder et al. 
2001) and two MOS CCD (Turner et al. 2001) cameras as 
the imaging spectrometers. All XMM-Newton instruments 
operate simultaneously. XMM-Newton observed MCG-5- 
23-16 twice on 13 May 2001 and 1 December 2001 for 38ks 
and 25ks, respectively. The EPIC PN and MOS observa- 
tions were carried out in the full frame mode using the 
medium filter. Here we present the data from the PN and 
MOS cameras. 

The raw PN and MOS events were processed and fil- 
tered using the most recent updated calibration database 
and analysis software (SAS v5.3.3) available in Decem- 
ber 2002. Events in the bad pixels file and those adja- 
cent pixels were discarded. Events with pattern 0-4 (sin- 
gle and double) for the PN and 0-12 (similar to ASCA 
event grades 0-4) for the MOS cameras were selected for 
both the observations of 2001. Examination of background 
light curves extracted from source-free regions showed that 
both the observations were affected by particle induced 
flares. These flares are characterized by strong and rapid 
variability in the background light curve. The particle 
induced events were filtered out by excluding the time pe- 
riods where the count rate increased by 3a from the quies- 



cent state background rate. This resulted in net exposure 
times of 8.5ks and 13. Iks for the PN and MOS, respec- 
tively for the observation of May 2001. The observation of 
December 2001 was not affected severely by the particle 
induced flares: the net exposure times are 19.7ks for the 
PN and MOS cameras. 

The source spectra were extracted from the final filtered 
event lists using a circular region of radius 90" centered on 
the observed position of MCG-5-23-16 for both the PN and 
MOS cameras. Background spectra were extracted using 
appropriate annular regions around the observed position 
of MCG-5-23-16. Any source in the background region was 
masked out. Appropriate response and effective area files 
for both the PN and MOS cameras were created using the 
SAS. 

3. SPECTRAL ANALYSIS 

The PN and MOS spectra of MCG-5-23-16 were ana- 
lyzed using the spectral fitting package XSPEC vl 1 . 2 . 0. In 
order to check for consistency in the spectral calibration of 
PN and MOS cameras, to check for spectral variability and 
to derive time averaged spectral characteristics of MCG- 
5-23-16, we have carried out spectral model fitting (i) to 
the individual PN and MOS data sets extracted from each 
observation, (ii) jointly to the PN and MOS data sets for 
each observation, and (iii) to the time averaged PN data 
obtained by combining the two PN data sets from the two 
observations. Thus, for spectral fitting, we formed nine 
data sets from the two observations carried out using the 
three imaging cameras: PN, MOS1, and MOS2. The nine 
data sets were first fitted by a simple absorbed power-law 
model over the entire 0.2 — 12 keV bandpass. None of the 
fits was acceptable e.g., the power-law fit to the combined 
PN spectrum resulted in a minimum \ 2 {Xmin) °f 2503.7 
for 504 degrees of freedom (dof). To show the significant 
deviations from a best-fit power law and to determine the 
true continuum shape, we refitted the absorbed power-law 
model after excluding the significant features. For these 
fits we used the data in the 2.5 — 5 keV and 7.5 — 10 keV 
bands. In these bands, AGNs usually do not show any 
significant localized feature. Figure 1 shows the combined 
PN data fitted with the absorbed power-law model and the 
ratio of data to the best-fit model extrapolated to lower 
and higher energies. As can be seen in Fig. 1, the poor 
fit is due to the significant features: a soft excess com- 
ponent below ~ 1 keV, a strong iron Ka emission line at 
~ 6.4 keV, and likely excess emission above 10 keV. This 
paper concentrates on the 2.5 — 12 keV band and the soft 
X-ray spectrum will be studied in a separate paper using 
RGS data also. 

3.1. The 2.5 — 10 kcV continuum 

To characterize the X-ray continuum, simple red-shifted 
power-law models modified by absorption due to an inter- 
vening medium at z = were fitted to the six spectra 
derived from the two observations with the three EPIC 
cameras: PN, MOS1, and MOS2. These fits were carried 
out over the 2.5 — 10 keV band excluding the FeKa re- 
gion (5 — 7.5 kcV). All the fits resulted in acceptable % 2 ; 
however, the best-fit values of equivalent hydrogen column 
(A%) were much higher (Nr ~ 1.7 x 10 22 cm" 2 for the 
PN data of December 2001) than the Galactic column of 
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8.82 x 10 20 cmT 2 (Elvis, Wilkes, & Lockman 1989). Set- 
ting the iVn parameter to the Galactic value resulted in 
unacceptable fits for all the data sets, thus implying heavy 
obscuration. The excess absorption column was then de- 
termined by introducing an additional absorption compo- 
nent at the source redshift and by refitting the power-law 
model to individual data sets. Table 1 lists the best-fit 
parameters: excess Nr, the photon index Fx, and the fit- 
statistic Xmin- AH t ne errors quoted, here and below, are 
at the 90% confidence level. Also listed in Table 1 are the 
observed flux in the 2.5 — 10 keV band and the absorp- 
tion corrected (intrinsic) flux in the same band. As can be 
seen in Table 1, the best-fit parameters Nnand Fx, derived 
from the data obtained simultaneously with different in- 
struments, are consistent within the errors. However, the 
observed flux is found to differ by ~ 20% between EPIC 
PN and MOS cameras. This discrepancy is not unusual 
for X-ray instruments in their early phase of calibration. 

In order to further constrain the spectral shape, the ab- 
sorbed power-law model was fitted to the PN and MOS 
spectra jointly. The relative normalizations for the dif- 
ferent instruments were kept free allowing for the small 
difference in the calibrated absolute flux, and any differ- 
ences in the fraction of encircled counts contained in the 
PN and MOS extraction cells. The results of these fits 
are also listed in Table 1. There are small variations in 
the best-fit spectral parameters between the two observa- 
tions: AN K = (4.4±2.7) x 10 21 cm" 2 ; Ar x = 0.13±0.06; 
A/j„t = 7 x 10~ 12 erg s _1 cm~ 2 . The source became 
harder and fainter in Dec. 2001. 

To determine the average continuum shape, we fitted 
the absorbed and red-shifted power-law model to the time 
averaged PN data obtained by combining the two PN data 
sets from the two observations. The time averaged photon 
index is Tx = 1-69 ± 0.03 and the observed flux in the 
2.5 - 10 keV band is 7 x 10~ n erg s" 1 cm" 2 . 

3.2. The iron Ka emission 

Figure 2 shows the ratio of the PN spectrum extracted 
from the December 2001 observation to the best-fit power 
law derived by using the data in the 2.5 — 5 keV and 
7.5 — 10 keV bands as described in the preceding section 
(see Table 1). This plot clearly reveals a broad emission 
feature in the 5 — 7 keV band and peaking at ~ 6.4 keV. 
This feature is common among Seyfert galaxies and is at- 
tributed to the fluorescence emission from the iron K-shcll 
(e.g., Nandra et al. 1997). The FeKa emission of MCG- 
5-23-16 appears to be nearly symmetric in its red and 
blue wings and is clearly different in shape from the red- 
shifted and asymmetric FeKa profiles generally observed 
from Seyfert 1 galaxies. Fig. 2 also suggests the likely 
presence of an iron K-edge at ~ 7.1 keV. To measure the 
strength and shape of the FeKa emission, we parameterize 
the observed profile in terms of simple Gaussian models. 
For this analysis, the co-added PN data were used, as the 
PN has the largest effective area among the three EPIC 
cameras. Hence the derived parameters characterize the 
time averaged FeKa profile. Initially, a narrow unresolved 
Gaussian (a = 0.01 keV) was used for the FeKa emission 
and the absorbed and red-shifted power-law model was 
used for the continuum as described in section 3.1 (Ta- 
ble 1). This model resulted in a poor fit with Xmin 01 



262 for 235 dof. The ratio of the data to the best-fit nar- 
row Gaussian model is plotted in Figure 3 (second panel 
from the top) which shows a dip at ~ 6.4 keV and excess 
counts on the red and blue sides of the dip. This suggests 
a strong core at ~ 6.4 keV and broad wings on both sides 
of the core. Varying the width of the Gaussian improved 
the fit significantly (A% 2 = 49 for 1 additional parameter). 
This suggests that the line may be broad. Addition of a 
narrow Gaussian component (a = 0.01 keV) at ~ 6.4 keV 
improves the fit significantly ((A^ 2 = 26.3 for two addi- 
tional parameters). The best-fit parameters describing the 
observed FeKa profile are listed in Table 2. The peak en- 
ergies of both the narrow and broad components are con- 
sistent with neutral iron. The full width at half maximum 
(FWHM) of the broad component is ~ 42000 km s" 1 . Ad- 
dition of an absorption edge at 7. 1 keV results in marginal 
improvement in the fit (A% 2 = 2.7 for one additional pa- 
rameter) at a significance level of <~ 93% based on an F-test 
(Bevington 1969). 

The width of the broad iron Ka line observed from 
MCG-5-23-16 is too large to be produced in regions other 
than the accretion disk around a super-massive black hole 
(SMBH). Therefore we checked whether the observed pro- 
file is consistent with a relativistic disk-line model (Fabian 
et al. 1989). This model assumes a Schwarzschild geome- 
try and the disk emissivity is a power-law function of disk 
radius r i.e. k r'. We fixed the inner disk radius at 6r g and 
the outer disk radius at 500r g , where r g — GM/c 2 . First 
we fitted the disk-line model without a narrow Gaussian 
component. The continuum was the absorbed power-law 
model as before. The free parameters were the disk in- 
clination angle (i) between our line of sight and the disk 
normal, the disk emissivity q, and the normalization of 
the disk line. This fit resulted in an acceptable fit statistic 
{ Xmin = 229.8 for 236 dof). The ratio of the data and the 
best-fit model is plotted in Figure f4 (top panel). The dip 
seen at ~ 6.4 keV is due to the characteristic shape of the 
line profile from a face-on disk - namely a strong core and 
a red wing. The disk inclination angle was found to be 
in the range 0° — 11°. This is physically inconsistent with 
a Seyfert 1.9 galaxy according to the standard unification 
scheme of Seyfert galaxies (Antonucci 1993). In view of 
the fact that the observed FeKa profile can be described 
by a combination of narrow and broad Gaussian and a 
narrow component is expected due to reflection from cold 
matter away from the disk, we added a narrow Gaussian 
component to our disk-line model. The energy of the nar- 
row component was fixed at 6.4 keV. This fit resulted in 
significant improvement over that without a narrow com- 
ponent (Ax 2 = 32.4 for one additional parameter). This 
is an improvement at a significance level of > 99.99%. The 
ratio of the data to the best-fit model is shown in Figure 4 
(second panel from the top). The best-fit inclination angle 
now is i = 46.3^38- The equivalent width of the disk-line 

is 129^23 e V and the line energy is consistent with neutral 
iron. Addition of an absorption edge at 7.1 keV does not 
lead to a significant improvement in the fit, and the upper 
limit to the optical depth is found to be r < 0.03. 

We have also tested the observed iron Ka emission pro- 
file of MCG-5-23-16 against that expected from an accre- 
tion disk around a rotating (Kerr) black hole (Laor 1991). 
We fixed the inner radius at 1.23r s , the minimum ra- 
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dius allowed for a Kerr black hole, and outer radius at 
400r 3 which is the maximum radius allowed by the model. 
As before, the disk emissivity is a power-law in radius 
(oc r~ q ). The variable parameters were EpeKa, h and 
the line normalization. We used an absorbed and red- 
shifted power-law model for the continuum. The bottom 
two panels in Fig. 4 show the ratio of the data and best-fit 
disk-line models for a Kerr black hole. Without a narrow 
Gaussian component, the fit is acceptable (x 2 = 251.7 for 
235) dof ) but the best-fit inclination angle is ~ 18 deg. In- 
clusion of a narrow Gaussian component (a = 0.01 keV) 
resulted in a significant improvement (Ax 2 = 60.5 for one 
additional parameter) and in the best-fit inclination an- 
gle of ~ 47.4 deg (see Table 2). Addition of an absorp- 
tion edge at 7.1 keV does not improve the fit significantly 
(A% 2 = 0.9 for one additional dof). 

In Figure 5, we have plotted the final best-fit models 
for the observed iron Ka line profile of MCG-5-23-16. The 
three models - (i) a combination of a narrow Gaussian and 
broad Gaussian, (ii) a disk line for a Schwarzschild or (iii) 
a Kerr black hole all describe the observed data equally 
well. 

In order to look for any variations in the strength of 
the FeKa emission, we have also analyzed the observed 
FeKa profiles from the two XMM-Newton observations. 
We used both PN and MOS data for this purpose. We 
fitted a combination of a broad and a narrow Gaussian 
profile to the PN and MOS data jointly for each observa- 
tion. The 2.5 — 10 keV continuum was fixed as determined 
earlier (see Table 1). The rest-frame line energies of both 
the narrow and broad components were fixed at 6.4 keV. 
The results of these fits are listed in Table 3 for both the 
observations. It appears that the strength of the broad 
component of the FeKa line decreased by about a factor of 
two between the two observations taken six months apart, 
while the narrow component did not vary significantly. 

3.3. The Compton reflection model 

Under certain conditions, the Compton reflection emis- 
sion can be important even below ~ 10 keV (see sec- 
tion 4.2). If the onset of the reflection component is 
at 5 keV, the presence of a strong Fe K-shell edge at 
~ 7.1 keV may artificially mimic a broad emission line fea- 
ture at ~ 6.4 keV. This artificial "broad line" is likely to 
arise if the observed continuum is flatter and is dominated 
by the reflection component. The 2.5 — 10 keV photon in- 
dex of MCG-5-23-16 derived from the combined PN data 
(Tx ~ 1-7) is flatter than the average 2 — 10 keV photon 
index of Seyfert galaxies (< T x >~ 1.9; Mushotzky 1997). 
This flatter spectrum and strong iron Ka line of MCG-5- 
23-16 could then be partly due to the contribution of a 
Compton reflection process in a geometry where the pri- 
mary X-ray source is partially covered by cold matter e.g. 
a torus. The reflection emission has been detected from 
MCG-5-23-16 with RXTE observations in 1995 (Weaver, 
Krolik, & Pier 1998). To test whether the broad fea- 
ture near 6.4 keV could be due to reflection and iron K- 
edge absorption of the primary emission, we have fitted a 
Compton reflection model to the combined PN data in the 
2.5 — 12 keV band. We used the reflection model "pexrav" 
(Magdziarz & Zdziarski 1995) available with XSPEC. This 
model calculates the reflected spectrum from a neutral disk 



exposed to an exponentially cut-off power-law spectrum 
(Magdziarz & Zdziarski 1995). Since XMM-Newton does 
not cover the expected peak (~ 30 — 40 keV) of the reflec- 
tion component, it is not possible to constrain all the pa- 
rameters of the reflection model using XMM-Newton data 
alone. Instead we fix the cut-off energy of the primary 
power-law at 200 keV, disk inclination at 50 deg, and the 
abundance of heavy elements at the solar value. The free 
parameters are the photon index of the primary power-law, 
the power-law normalization, and the relative amount of 
reflection compared to the directly viewed primary spec- 
trum (R). We also included a narrow Gaussian line model 
to describe the FeKa emission. The results of this fit are 
listed in Table 4. The best-fit photon index of the incident 
power-law continuum is ~ 1.8 and is consistent with the 
mean photon index for Seyfert 1 galaxies. The ratio of 
the observed data to the best-fit model is shown in Fig- 
ure 6. Although the fit is acceptable (Xmin = 224.4 for 
269 dof), there is weak excess emission near 6 keV sugges- 
tive of the possible presence of a weak broad emission line. 
To investigate further, we included an additional Gaus- 
sian component in our reflection model and carried out 
the fitting. This resulted in an improvement in the fit 
(Ax 2 = 20.1 at the cost of 2 dof). This is an improvement 
at a significance level of 99% according to the maximum 
likelihood test. The best-fit parameters are listed in Ta- 
ble 4. The relative normalization of the reflection emission 
derived from the PN data is much higher than the value of 
R = 0.54±°;f 7 inferred from the BeppoSAX observations 
of MCG-5-23-16 (Risaliti 2002). This is due to the excess 
emission above ~ 10 keV which amounts <~ 20% above 
10 keV (see Fig. 1). Such a large discrepancy in the re- 
flection emission inferred from the PN and Bepposax data 
could be due to calibration uncertainties. This view is sup- 
ported by MOS data, albeit lower signal-to-noise, which do 
not show a clear excess emission above 10 keV. 

To test the reliability of the excess emission above 
10 keV in the XMM-Newton PN data, we have analyzed 
the EPIC PN spectrum of a well known BL Lac, Mkn 421. 
The BL Lac type AGNs do not show the Compton re- 
flection component, just a smooth continuum. Therefore, 
they are good targets to check the reliability of spectral 
features in other AGN X-ray spectra. We retrieved the 
data corresponding to the XMM-Newton observation of 
Mkn 421 on 4 May 2002 for an exposure time of ~ 40 ks 
from the public archive maintained at HE AS ARC. The 
data were processed in a similar way as described above for 
MCG-5-23-16. The source spectrum was extracted from 
an annular region with inner and outer radii of 10" and 
80", respectively and centered at the observed source posi- 
tion. Due to the high flux of Mkn 421, the events from the 
core of the point spread function were excluded to avoid 
possible pile-up which results in artificial hardening of the 
spectrum. The background spectrum was extracted from 
the source free regions near the position of Mkn 421. New 
response files were created for Mkn 421. The PN spectrum 
was then grouped appropriately and fitted by an absorbed 
power law model in the 2.5 — 10 keV band. Figure 7 shows 
the PN spectrum of Mkn 421 fitted with a power-law and 
the ratio of data and the best-fit power-law model. An 
excess over the best-fit power-law is seen above 10 keV. 
This excess is ~ 20 — 30% above ~ 10 keV which is similar 
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to that seen in the PN spectrum of MCG-5-23-16 (Fig. 1). 
Since Mkn 421, being a blazar, is not expected to show 
any excess emission, we conclude that the excess emission 
over the power-law seen above ~ 10 keV in the spectra of 
Mkn 421 and MCG-5-23-16 is likely an artifact and reflects 
poor spectral calibration of the PN above 10 keV. Lack 
of any significant feature in the 2.5 — 10 keV band of the 
Mkn 421 PN spectrum suggests that the spectral calibra- 
tion of PN in that energy band is indeed good and the 
presence of a broad iron Ka emission is not an artifact. 

We refitted the Compton reflection model to the com- 
bined PN data after excluding the data above 10 keV. We 
utilized the best-fit parameters for the reflection compo- 
nent derived from the BeppoSAX observations by Risaliti 
(2002) and fixed the cut-off energy Ec at 157 keV and rel- 
ative normalization R at 0.54. The best-fit parameters are 
listed in Table 4. It is clear that a narrow Gaussian and 
reflection emission cannot account for the strength of the 
broad feature at ~ 6.4 keV which implies the presence of 
a broad FeKa line. 

4. DISCUSSION 

We have presented the results from the two observations 
of MCG-5-23-16 using XMM-Newton. The 2-10 keV flux 
is found to be 6.5 x 10~ n erg s _1 cm~ 2 during December 
2001 and 7 x 10~ n erg s _1 cm~ 2 during the June 2001 ob- 
servation, comparable with the 7.3 x 10~ n erg s _1 cm~ 2 
measured during the ASCA observations of 1994 (Weaver 
et al. 1997). The 2.5 — 10 keV photon index appears 
to have flattened slightly from 1.77±°:o5 ( Ma y 200 1 ) to 
1-641q q3 (December 2001). Both these photon indices 
are steeper than the index of 1.4 — 1.5 observed with 
Gmgawhen the source was in a lower flux state (Nandra & 
Pounds 1994; Smith & Done 1996). The X-ray spectrum of 
December 2001 is slight flatter than the ASCA spectrum of 
1994 when the source flux was slightly higher. These obser- 
vations suggest that the X-ray spectrum flattens at lower 
flux levels - a general trend well established for Seyfert 
galaxies (Singh, Rao, & Vahia 1991; Done, Madejski, & 
Zdotycki 2000; Zdziarski & Grandi 2001; Vaughan & Edel- 
son 2001; Nandra 2001; Dewangan et al. 2002). 

4.1. FeKa emission and its variability 

XMM-Newton observations of MCG-5-23-16 have reaf- 
firmed the presence of a complex iron Ka emission ear- 
lier observed with ASCA (Weaver et al. 1997; Weaver & 
Reynolds 1998; Turner, George, Nandra, & Mushotzky 
1998) and RXTE (Weaver, Krolik, & Pier 1998). Our spec- 
tral fits imply that the line is broad, superimposed with 
a narrow unresolved component, and possibly affected by 
an absorption edge of iron at 7.1 keV. The width of the 
broad line derived by fitting a Gaussian profile suggests 
that this component could only be formed in an accretion 
disk around a black hole. Indeed the broad component 
is consistent with that expected from a relativistic accre- 
tion disk around a Schwarzschild or a Kerr black hole. 
Our results also confirm the ASCA result that it is dif- 
ficult to explain the entire line profile as originating en- 
tirely in an accretion disk. The Gaussian and relativistic 
line model fits to the iron Ka line imply two distinct com- 
ponents. The narrow component is unresolved (FWHM 
< 1000 km s _1 ) and has an equivalent width of ~ 40 cV. 



Similar narrow unresolved components have been observed 
with Chandra and XMM-Newton from a number of AGNs 
(e.g., Kaspi et al. 2001; Pounds et al. 2001; Turner et al. 
2002), and are found to have their peak energy at 6.4 keV 
with equivalent widths of ~ 50 eV . Such narrow com- 
ponents are naturally expected from the torus and/or the 
broad line region in the framework of the unification model 
of Seyfert galaxies. The maximum equivalent width of an 
FeKa line produced by fluorescence in a cold reflector far 
from the disk with iVn — 10 23 cm~ 2 is ~ 65 eV in the case 
of NGC 4258 (Reynolds, Nowak, & Maloney 2000). Thus 
the observed equivalent width of the narrow iron Fa line 
of MCG-5-23-16 is consistent with that expected from the 
torus. 

The broad component has an FWHM of ~ 4 x 
10 4 km s" 1 and an equivalent width of ~ 110 eV. This 
component is described by a broad Gaussian or rela- 
tivistic disk-line arising from an accretion disk around a 
Schwarzschild or Kerr black hole suggesting that the line 
profile is not strongly modified by gravitational effects but 
broadened by the Doppler effect in a rotating disk. The 
disk-line fits imply that the disk is neutral and the disk- 
normal is inclined at an angle of ~ 50 deg. The likely 
presence of an absorption edge at ~ 7.1 keV (see Fig. 2) 
also suggests that the disk is neutral. The best-fit values 
of the emissivity indices inferred from the disk-line fits (in 
Schwarzschild and Kerr geometry) are similar but lower 
than the average value of 2.5 obtained for Seyfert 1 galax- 
ies (Nandra et al. 1997). It should be noted that Nandra 
et al. had also fixed the inner radii to 6r g (Schwarzschild 
black hole) and 1.23r ff (Kerr black hole). Weaver et al. 
(1997) inferred q to be in the range of 5 — 10 but allowed 
ri to vary with the best-fit values in the range 9 — 22. 
These results imply that if the emissivity law for MCG-5- 
23-16 is the same as that of an average Seyfert 1 galaxy, 
then the inner radius of the accretion disk is larger for 
MCG-5-23-16. 

Our spectral fits suggest that the 2.5—10 keV continuum 
has flattened from T x ~ 1.77 (May 2001) to T x ~ 1.67 
(December 2001) and the broad component of the iron Ka 
line has weakened from an equivalent width of ~ 180 eV 
(May 2001) to - 90 cV (December 2001). In the ASCA 
observations of May 1991, the photon index and the equiv- 
alent width of the broad component of the iron Ka line 
were found to be 1.79 ± 0.08 and 254+yy, respectively. 
These results suggest weakening of the broad iron Ka line 
with flattening of the X-ray continuum. Higher equiva- 
lent widths of the iron Ka emission could be produced by 
an increased abundance of iron or by increasing the re- 
flection component. It is unlikely that the change in the 
equivalent width of iron Ka is due to a change in the iron 
abundance due to the following reasons: (i) it is very un- 
likely that iron abundance changes by a factor of about 
two in just six months, and (ii) it is difficult to explain 
the variation in the shape of the non-thermal X-ray con- 
tinuum by changing the iron abundance. Increasing the 
reflection emission from a neutral disk may flatten the 
continuum. This will require a change in the geometry 
of the accretion disk-corona system (e.g. increasing the 
covering factor at the corona will result in more reflec- 
tion emission). However, this will increase the equivalent 
width of the iron Ka line with flattening of the contin- 
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uum. Large equivalent widths of FeKa can be produced 
if the accretion disk is ionized (Weaver et al. 1997). De- 
wangan (2002) reported a correlation between the shape 
of the X-ray spectrum and the energy of the FeKa line 
among radio-quiet type 1 AGNs and suggested that the 
accretion disks of AGNs with steeper continua are more 
ionized. The trend in Fx and the equivalent width of iron 
Ka emission from MCG-5-23-16 could be understood if the 
shape of the X-ray continuum is coupled with the ioniza- 
tion stage of the disk material in MCG-5-23-16 in a similar 
way to that observed among the radio-quiet AGNs. 

4.2. Is the broad feature at <~ 6.4 keV really an iron Ka 
line from an accretion disk? 

As shown in section 3, the broad feature at 6.4 keV can 
be modeled broad Gaussian (FWHM~ 40000 km s _1 ) 
or a line profile arising from an accretion disk around a 
Schwarzschild or a Kerr black hole. Alternatively, most of 
the flux in the broad component can also be accounted by 
reflection emission . However, the derived reflection frac- 
tion (R ~ 3) is much larger than that inferred from the 
BeppoSAX observations of MCG-5-23-16 (R ~ 0.5; Risal- 
iti (2002)) and the average value for Seyfert 1 galaxies 
(R <~ 1; Perola et al. 2002). This is not entirely impos- 
sible if the reflector is a putative torus far away from the 
central continuum source. The reflection fraction can be 
large if the central source was brighter before a time pe- 
riod which is equal to the time delay between the reflected 
and primary X-ray emission. However, such a large re- 
flection fraction for MCG-5-23-16 is unlikely due to the 
following reasons - (i) the calibration of the EPIC PN 
above 10 keV is reliable at a level of ~ 20% (Snowden, 
private communication). MOS data do not show clear ex- 
cess emission above 10 keV as is the case for PN. (ii) The 
broad feature at ~ 6.4 keV weakens with flattening of the 
continuum. This is not expected if the broad feature were 
the reflection emission. Simultaneous observations with 
XMM-Newton and Integral would be crucial to character- 
ize both the FeKa emission and reflection component. 

4.3. The Iron K-shell Edge 

XMM-Newton observations of December 2001 have re- 
vealed the presence of an Iron K-shell absorption edge near 
~ 7.1 keV. However, the edge is not well constrained and 
is detected only at a significance level of 93%. The edge en- 
ergy is not well determined but it appears to be ~ 7.1 keV 
(see Fig. 2). The optical depth is found to be 0.03±°;°^. 
Assuming that the edge arises due to the K-shell absorp- 
tion of neutral iron and the corresponding photoionization 
cross section (2.6 x 10~ 20 cm~ 2 ; Donnelly, Bell, Scott, & 
Keenan 2000), we determine the column density of neutral 
iron to be Np c « 1.1 x 10 18 cm~ 2 . This corresponds to an 
effective Hydrogen column density of Ah ~ 4 x 10 22 cm~ 2 , 
which is a factor of two higher than the hydrogen column 
density derived by fitting an absorbed power-law model to 
the observed data. Given the uncertainty in the derived 
value of the optical depth, it cannot be ruled out that the 
Fe K-edge arises entirely in the absorbing material along 
our line of sight in the form of a putative torus. It is inter- 
esting to investigate whether both the absorption K-edge 
and the narrow component of the Ka line due to iron could 
be produced in the the putative torus around an accretion 



disk. In the case of isotropic emission, the equivalent width 
of the emission line should be ~ QAfcN^Ape x 10~ 23 . For 
a typical covering factor f c ~ 0.7 for the torus and solar 
abundances (Af e = 1), the absorption column required to 
produce an equivalent width of 40 eV for the narrow iron 
Ka line is ~ 10 23 cm~ 2 . This column is larger by a fac- 
tor of five than the absorption column derived from the 
power-law fit. It is also larger by about a factor of three 
than the absorption column derived from the iron K-edge, 
since the absorption features are produced only by a col- 
umn along our line of sight while the emission features are 
spatially integrated in the case of a source which is unre- 
solved. The above results imply that we are viewing the 
outer edge of a putative torus where the column density is 
lower than the inner regions of the torus. This is consis- 
tent with that expected for a Seyfert 1.9 galaxy from the 
unification scheme (Antonucci 1993). 

5. CONCLUSIONS 

We studied the 2.5 - 10 keV spectrum of MCG-5-23-16 
which was observed twice with XMM-Newton in May and 
December 2001. Our main results are as follows. 

1. The iron Ka line profile of MCG-5-23-16 is broad 
and nearly symmetric. The observed profile 
could be modeled by a combination of a nar- 
row (unresolved) and a broad Gaussian compo- 
nent. The broad component is also consistent with 
that expected from an accretion disk around a 
Schwarzschild or a Kerr black hole. 

2. The disk inclination angle of ~ 46°, inferred from 
the disk- line fits, is consistent with that expected 
for a Seyfert 1.9 galaxy according to the Seyfert 
unification picture. 

3. Alternatively, most of the flux of the broad feature 
at 6.4 keV can also be modeled as the reflection 
emission, however, the reflection fraction (~ 3) is 
much larger than that derived from the BeppoSAX 
observations (R <~ 0.5). This is likely to be due to 
uncertainty in the calibration of EPIC PN above 
10 keV. 

4. We have derived three values of absorption col- 
umn (i) N-£i ~ 10 22 cm~ 2 from the best-fit ab- 
sorbed and red shifted power-law model, (ii) ~ 
4 x 10 22 cm~ 2 from the absorption edge, and (Hi) 
Ah ~ 10 23 cm~ 2 required to produce the narrow 
iron Ka line emission. The above three values are 
consistent with a geometrical picture in which our 
line of sight passes through the outer edges of a pu- 
tative cold torus obscuring the nuclear X-ray source 
as expected for a Seyfert 1.9 galaxy from the unifi- 
cation scheme. 

5. The broad iron Ka emission appears to weaken with 
the flattening of the X-ray continuum. This can be 
understood if the degree of ionization of iron in the 
accretion disk decreases with flattening in the X- 
ray spectrum. However, in any case the disk does 
not appear to be highly ionized (i.e. Hydrogen or 
Helium- like iron). 
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Fig. 1. — Ratio of combined PN data and the best-fit model derived by fitting an absorbed and red-shifted power-law to the observed data 
over the energy band of 2.5 — 10 keV but excluding the iron K« region 5.0 — 7.5 keV band. A strong soft X-ray excess below 2 keV and 
an emission feature at ~ 6.4 keV due to iron K-shcll emission are seen clearly. Also seen is the excess emission over the power law above 
~ 10 keV. 
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Fig. 2. — The FcKa region of the X-ray spectrum of MCG-5-23-16 shown as the ratio of the observed data of December 2001 and the 
best-fit power-law model. The best-fit power-law continuum was derived by using the data in the 2.5 — 5 keV and 7.5 — 10 keV bands. A 
strong and broad FcKo peaking at ~ 6.4 keV is evident. 
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Fig. 3. — Ratios of combined PN data to the best-fit Gaussian models of the iron Ka line profile. From top, the best-fit models are ( 
absorbed and red-shifted power law (PL), (ii) PL and a narrow Gaussian with a = O.Of keV, (iii) PL and a combination of a narrow and 
broad Gaussian, (iv) same as above but with an absorption edge at 7.1 keV. 
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Fig. 4. — Ratios of combined PN data to the best-fit relativistic disk-line models for the iron K« line profile. From the top, the best-fit 
models are (i) an absorbed and red-shifted power-law and a disk-line model for a Schwarzschild black hole, (ii) same as above with an addition 
of a narrow Gaussian line with a = O.Of keV, (iii) PL and disk-line model for a Kerr black hole, (iv) same as above with the addition of a 
narrow Gaussian (a = O.Of keV). 
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Fig. 7. — The EPIC PN spectrum of Mkn 421 and the best-fit power-law model over the 2.5 — 12 kcV band. A clear excess over the best-fit 
power-law is seen at energies above ~ 10 keV. This excess is similar within error to that seen in the PN spectrum of MCG-5-23-16. 
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Table 1 

Best-fit absorbed power-law model parameters for MCG-5-23-16 



Observation 


Data a 


Model b 


N H C 


Tx 


fobs 


font 


xVdof 








10 22 cm" 2 




2001 December 


PN 


PL 


1 67 +ui9 

1 -°'-0.19 


1 66 +,UB 

i,uu -0.03 


6.9 


7.6 


138.4/146 




MOS1 


PL 


i 9 o+0.42 
± -^°-0.42 


1 C9+0.09 
L -°*-0.09 


5.7 


6.1 


169.8/161 




MOS2 


PL 


i-56i°; 4 4 I 


1 cc+0.10 


5.5 


6.0 


166.2/161 




PN+MOS1+MOS2 


PL 


i-65t°;i! 




6.0 


6.5 


478.0/472 


2001 May 


PN 


PL 


1 OQ+0-28 
l.»J_ 28 


+0.04 
' d -0.04 


7.1 


7.9 


160.4/146 




MOS1 


PL 


9 n o+0.48 
z.uo_ 4g 


i-69±8;ig 


6.4 


7.2 


182.9/161 




MOS2 


PL 


1 89 +0 - 48 


1.70lS;ll 


6.2 


6.9 


195.7/161 




PN+MOS1+MOS2 


PL 


9 0Q+ 0,22 

z.u»_q 2i 


i 77 +0.04 


6.4 


7.2 


553/472 


2001 Dec. + May 


PN 


PL 


'"-0.17 


i fiq+ 03 


7.0 


7.7 


123.6/158 



"The data in the 2.5 — 10 keV were used. 

b Simple red-shifted power-law modified by absorption 

c Exccss absorption intrinsic to the source apart from that due to the Galactic column of 8.82 X 10 -20 cm -2 . 



Table 2 

Best-fit parameters of the two component Gaussian model for the iron Ka emission of MCG-5-23-16 



Parameter PL 

PN+MOS PN+MOS Combined PN 

Dec. 2001 May 2001 May and Dec. 2001 





NO EDGE 


EDGE 


NO EDGE 


EDGE 


NO EDGE 


EDGE 


r x 


1.64 


1.64 


1.77 


1.77 


±-O9_ .04 


1 65+ ,UJb 
1 - DO -0.06 




1.65 


1.65 


2.09 


2.09 


L - ,o -0.16 


1 60+ ' 22 
l.DU_ 22 


En 


6.4 


6.4 


6.4 


6.4 


6.4 


fi oq+0.02 
D-<J»_0.03 




0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


EW N 


48.0tl 8 2 


r- 9 0+6.6 
oz.o_ 14 5 


AO 1+13.2 

^■-'--le.o 


46 


™ o+9.6 
oo.o_ g g 


™ 9+9.6 


E B 


6.4 


6.4 


6.4 


6.4 


6 S2+ ' 07 


R 09+0. 09 
"■ oz -0.09 


OB 


o 7f) +0.109 


S76+ 0143 
u -°'°— 0.074 


0.550lS;iS 


0.602l°;i 8 | 


n Qoi +0.104 

U.OOly Q g7 


n orr+0.128 

U.OOO_Q QgQ 


EW B 


Q1 <;+21-5 
91.0_ 12 9 


92.6ti|1 


1 77+38 
L "-30 


189tf 7 


114+21.0 

ii4 -21.4 


96.4l 3 ; 6 7 


E e dge 




7.1 




7.1 




7.1 


T 




02+ - 02 
u - uz -0.02p 




02+ 03 
u ' uz -0.02p 




n3 +004 


x 2 ■ 


727.2 


722.9 


783.8 


782.5 


186.7 


184.0 


dof 


697 


696 


703 


702 


232 


231 
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Table 3 

Results for relativistic disk-line fits to the combined PN data 



Parameter 


Schwarzschild 


geometry 


Kerr g 


;eometry 




NO EDGE 


EDGE 


NO EDGE 


EDGE 


1 A. 


1.68 ±0.02 


1.68 ±0.02 


1.68 ± 0.03 


1 65+^ 


l\Tr»T"m ci 1 i f7 fi 1" i r ~i n ^ 
1\U1 IlldllZjajLlUll 


0.021 ±0.001 0.021 ±0.001 


n 091 -i- n nm 




iV H (xl0 22 cm- 2 ) 


i 79 +0.12 
lm ' z -0.16 


1 7 9 +0.14 
11 ' z -0.17 


1.72+2-i? 

— 0.21 


1.62+2-11 

—0.17 


E N ( keV) 


6.4 


6.4 


6.4 


6.4 


<tn( eV) 


100 


100 


100 


100 


£Wat( cV) 


43.5±|:§ 


44.71™ 


38.0+11 


38.5^-3 


£7 D ( keV) 


6-3i±S:Sg 


6.3ol8:S 


6 29+ 09 

°' zy -0.12 


6.30l°;? 9 2 


q 


i 01 +0.40 
J-- 81 - 0.37 


1 81 + - 41 
i - 8i -0.38 


i 7C+0.38 


i 7 o+0. 36 
L - ' J -0.61 


n/r g 


6 


6 


1.23 


1.23 


r /r g 


400 


400 


400 


400 


i 


4«-3±i:t 


46.0+j- 5 


47.4+^ 


47.9+?? 


EW D { eV) 




128l 2 J 


llS+jt 


llfi+27 
iiD -25 


E edge { keV) 








7.1 


r 




O.OlOl^ 1 




ns+ 08 


X 2 /dof 


190.7/232 


190.4/231 


191.1/232 


188.7/231 



a Power-law normalization at 1 keV in units of photons keV 1 s 1 cm 2 . 



Table 4 

Best-fit parameters of a reflection model for the continuum and Gaussian model for the iron Ka emission of 

MCG-5-23-16 



Parameter 




PL 4 


- Reflection 






Reflection parameters varied 


Reflection parameters fixed a 




Single Gaussian 


Double Gaussian Single Gaussian 


Double Gaussian 


N K (xlO 22 cm" 2 ) 


L - ii -0.2 


1 7C+U.O6 
i -'°-0.08 


1 7 rr+0.08 
J -'°-0.15 


1 67 +U ' 08 

i ' ' -0.05 


T x 


1 or+0.03 
i - 8o -0.05 


1 70+O.O2 
i -' S -0.02 


1 fi8+ - 03 


1 RQ+0-05 
i - b8 -0.02 


Efoid 


200 keV 


200 keV 


157 keV 


157 keV 


R 


9 c+0.3 


1 7+0.2 
J--' -0.2 


0.54 


0.54 


i 


60° 


60° 


60° 


60° 


En 


fi 07+0.03 
-0.02 


6.4 


6.4 


6.4 


on{ eV) 


0.01 


0.01 


0.01 


0.01 


EW N ( eV) 


52.0+°;? 


^ -^-11.4 


62.0t 9 6 j 


34.5li?- 2 


Eb 




6.2ot8:3 






ctb 




o.islS:S 




31 +0 12 


EW b 




90 9+33.5 
zy.z_ 13 g 




83.9l 2 2 J 


A,m%ri 


224.4 


204.3 


234.6 


184.3 


dof 


269 


267 


236 


233 



"Relative reflection R and cut-off energy Eq were fixed at the values derived from the BeppoSylAbbservations (Risaliti 2002). 



